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Abstract. - It is shown that the BBGKY equations for a particle interacting with ideal gas imply 
exact relations between probability distribution of path of the particle, its derivatives in respect 
to the gas density and irreducible many-particle correlations of gas atoms with the path. These 
relations visualize that the correlations of any order always significantly contribute to evolution 
of the path distribution, so that the exact statistical mechanics theory does not reduce to the 
classical kinetics even in the low-density (or Boltzmann-Grad) limit. 
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Introduction. — In the seminal book [1] Bogoliubov 
formulated powerful tools for statistical theory of trans- 
port processes and noises in many-particle systems, in par- 
ticular, the exact hierarchy of evolution equations (now re- 
ferred to as the Bogolyubov-Born-Green-Kirkwood-Yvon, 
or BBGKY, hierarchy [2]) for s -particle (s = 1,2,...) 
distribution functions. But in practice the theory was can- 
celed out by common prejudice that it must confirm, at 
least for dilute gases, the Boltzmann's "molecular chaos 
hypothesis" , - the "StoBzahlansatz" [3] , - and thus reduce 
to the Boltzmann equation. This idea is so much attrac- 
tive that Bogoliubov introduced own hypotheses which al- 
low to truncate the BBGKY hierarchy at s = 2 and per- 
form the desired reduction [1]. In fact, however, until now 
neither Boltzmann's nor Bogoliubov's assumptions have 
a rigorous substantiation based on the BBGKY hierar- 
chy itselfQ, and hence the Boltzmann equation also stays 
without substantiation^ . 

Here, I suggest a short visual proof of that all this never 
will be substantiated, even for arbitrary dilute gas (or 



1 We take in mind spatially inhomogeneous gas evolutions, of 
course, since in respect to strictly homogeneous (translation invari- 
ant) dilute gas the StoBzahlansatz is undoubtedly true. 

2 The frequently mentioned Lanford theorem [4] about gas of 
hard spheres under the Boltzmann-Grad limit concerns the so called 
"hard-sphere BBGKY hierarchy" (about it see also e.g. [2,5]) which 
is not a true BBGKY hierarchy since represents interactions of the 
spheres by invented terms like the Boltzmann collision integrals (i.e. 
postulates what should be proved, if any). 

We here consider regular smooth interaction potentials only and 
thus the true BBGKY equations where the interaction is represented 
by the forces (potential gradients) like in the parent Liouville equa- 
tion. But our conclusions naturally extend to the limit of singular 
"hard-sphere" interaction. 



Boltzmann-Grad gas). At that, in order to simplify the 
proof and at once strengthen it, instead of the usual gas I 
consider a test particle in thermodynamically equilibrium 
ideal gas whose atoms interact with this particle only but 
not with each other (thus concentrating on situation least 
favorable for inter-particle correlations). 

Our first step will be representation of the BBGKY hi- 
erarchy of equations for this system in terms of irreducible 
s-particle correlations between the test particle and s — 1 
gas atoms ( s = 1 corresponds to probability distribu- 
tion of variables of the test particle). Then, we discover 
that these equations produce a set of formally exact re- 
lations between any s-particle correlation and derivative 
of the previous (s — l)-particle one in respect to the gas 
density. Next, make sure that all the correlations keeps 
non-zero even under the low-density (or Boltzmann-Grad) 
limit and, consequently, no truncation of the BBGKY hi- 
erarchy can be made without qualitative detriment to its 
solution. Finally, we touch upon important properties of 
the correlations. 

The BBGKY equations and cumulant distribu- 
tion functions. — We are interested in random walk 
R(i) of a test molecule (TM) which is placed into ther- 
modynamically equilibrium gas and starts at initial time 
moment t = from a specified position R(0) = Ro ■ 

Let P and M denote momentum and mass of TM, 
m, Vj and pj [j = 1,2, ...) are masses, coordinates and 
momenta of gas atoms, $(r) is potential of (short-range 
repulsive) interaction between any of them and TM, and 
n is gas density (mean concentration of atoms). At arbi- 
trary time t > , full statistical description of this system 
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is presented by the chain of (k + l)-particle distribution 
functions (k = 0,1,2,...): Fo(t, R, P| Ro ;n) which is 
normalized (to unit) density of probability distribution of 
TM's variables, and F k (t, R, r (fc ), P, p( fc) | R ;n) (where 
r ( fc ) — {rj... rj. } , p' fc ) ~ {pi... Pfc } ) which are probabil- 
ity densities of finding TM at point R with momentum 
P and simultaneously finding out some k atoms at points 
Tj with momenta p^ . A rigorous definition of such distri- 
bution functions (DF) was done in [1]. In respect to the 
coordinates r, they are not normalized, but instead (as 
in [1]) satisfy the conditions of vanishing of inter-particle 
correlations under spatial separation of particles (in other 
words, DF satisfy a cluster property with respect to spa- 
cial variables). Subject to the symmetry of DF in respect 
to Xj — {ijjPj} these conditions can be compactly writ- 
ten as follows: F k — > -Ffe-i G m (pk) at r k — > oo , where 
G m (jp) is the Maxwell momentum distribution of a parti- 
cle with mass m . 

The enumerated DF obey the standard equations [1]: 
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dt 
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r fc+i) Fk+i 
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with k — 0, 1, ... and along with obvious initial conditions 



Fk\t=o : 
6(R - Ro) 



<5(R 



exp(-H k /T) = 
-.1 - R) G m ( Pj ) 



(2) 



is Hamiltonian of subsystem " k atoms 



means the Poisson brackets, J, . 



where Hk 
TM", [..., 

fj...dr k dp k , $'(r) = V$(r) , and E(r) 
exp [— $(r)/T] . Notice that one can treat TM as molecule 
of non- uniformly distributed impurity and JT]) as equations 
of two-component gas [1] in the limit of infinitely rare im- 
purity with main component being in spatially uniform 
and thermodynamically equilibrium state. 

Equations ([I]) together with ((2|) unambiguously deter- 
mine evolution of Fq and eventually probability distri- 
bution of total TM's displacement, or path, R R . 
These equations will become more transparent if we make 
a proper linear change of DF Fk by new functions V k , 
namely, with the help of recurrent relations as follow: 



F (i,R,P|R ;n) = V (t, R, P| R ; ? 
fi(i,R,ri,P,p 1 |R ;n) = 
= V r (t,R,P|R ;n)/(ri-R,pi)- 
+ V r i(t,R,ri,P,p 1 |R ;n) , 

where /(r,p) = E(r)G m (p) , 

F 2 (t,R,rW,P,pW\Ro;n) -- 
= F (t,R,P|Ro;n)/(p 1 ,pi)/(p 2 ,p 2 ) - 

+ Fi(t,R,ri,P,pi|Ro;n)/(p 2 ,P2) - 
+ Pi(t,R,r 2 ,P,p 2 |R ;n)/(pi,pi) - 
+ ^ 2 (t,R,r( 2 ),P,p( 2 )|R ;n) 



(3) 



Apparently, from viewpoint of the probability theory, 
Vk are a kind of cumulants, or semi- invariants, or cu- 
mulant functions (CF). It is important to notice that if 
all these CF were zeros then all conditional DF of gas, 
Fk/Fo , would be independent on initial position R of 
TM and thus on its displacement R R . This fact makes 
visible very interesting speciality of the CF Vk : they are 
not mere correlations between instant dynamic states of 
TM and k gas atoms but simultaneously their irreducible 
correlations with the total past TM's displacement. 

Relation of many-particle correlations to prob- 
ability law of diffusion of the test particle. In 

terms of CF the BBGKY hierarchy acquires tridiagonal 
structure (we omit uninteresting algebraic details): 

BV Fi r 

' k = [H k , V k ]+n— $'(R - r k+1 )V k +i + 

°r Jk+l 



dt 
k 



TJ2 P kj G m (p k )E'(r k -R) 

3=1 
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MT dP 



Vk-i . (4) 



Here £'(r) = V£(r) , and Pkj symbolizes transposition 
of the pairs of arguments Xj and Xk ■ At that, initial 
conditions @ and the above-mentioned clustering condi- 
tions [1] take extremely simple form: 



V o (0 , R, P| R ; n) = S(R - R ) G M (P) , 
y fe (0,R,r( fe ),P,p( fe )|R o ;n) = , 
Vfc(i,R,r( fe ),P,pW|R ;n) -> at r^ ■, -> oo 



(5) 



( 1 < j < k). Thus, as it should be with cumulants, V k 
(k > 0) vanish under removal of already one of atoms. 

From equations ((4]) as combined with the boundary and 
initial conditions (JSJ) it is clear that passage to the limit 
in ([5]) must realize in an integrable way, so that integrals 
V k = J k+1 Vk+i take finite values. Let us consider them. 
By applying the operation J, to equations Q one obtains 



(6) 



+ T J2 F kjG m (p k ) E'{rk-R) 

3=1 
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V k - 



where 



Pi 



r 7 — R , and so on. 



(with k = 0,1,...). Because of ((5j) initial conditions to 
these equations are zero: V k {t = 0) = at any k . 

Now, in addition to V k , let us consider derivatives of CF 
in respect to the gas density, V k [ = dVk/dn . It is easy 
to see that differentiation of (0]) in respect to n yields 
equations for the V k which exactly coincide with ([6]) after 
changing there V k by V k . Besides, in view of ([5]), initial 
conditions to these equations again all are equal to zero: 
V k (t = 0) = at any k > . These observations strictly 
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imply the equalities Vj, = 14 , or 

^ F fe (i,R,r( fe ),P,p( fc )|R ;n) = (7) 

Jk+1 

This is main, original and formally exact, result of the 
present paper. Evidently, it confirms the assumed finitc- 
ness of integrals V k = J k+1 V k +i ■ Notice also that it an- 
ticipates similar but more complicated general statements 
of statistical kinetics of fluids. 

Discussion and resume. — The equalities (J7J con- 
tain the proof promised in the beginning of this paper. In- 
deed, they show, firstly, that all the many-particle correla- 
tions between gas atoms and total path, or displacement, 
of the test molecule (TM) really exist, i.e. differ from 
zero. Secondly, integral values of all the correlations, rep- 
resented in the natural dimensionless form, have roughly 
one and the same order of magnitude. Indeed, multiply- 
ing equalities by n and integrating them over TM's 
momentum and all gas variables, we have 

n k V k (t,A;n) = n k J...JJv k dP = 

k d k V (t,A;n) 
= n ^ CkVo(t,A;n) , 

where V (t,A;n) = J V (t, R, P|Ro; n) dP is just the 
probability density distribution of the TM's displacement, 
A = R — Ro , and c k some numeric coefficients obviously 
comparable with unit. Equivalently, unifying all CF into 
one generating function, we can write 

oc k k 

V {t, A;n) + Y, V k (t, A; n) = 
fc=i 

= V Q (t,A;(l + u)n) . (8) 

We see that distribution of total of the TM's random walk 
"is made of its correlations with gas atoms" like the walk 
itself is made of collisions with them. Hence, none of the 
correlations can be neglected if we aim at adequate anal- 
ysis of solution of the BBGKY equations. 

We see also that characteristic spacial volume occupied 
by the correlations has an order of the specific volume: 
(\J 1 ...J k JV k dP\/V ) 1/k ~ rT 1 . In the Boltzmann- 
Grad limit, n — > oo , ro — > (ro is effective radius of the 
interaction), 7rrQn = 1/A = const, this volume becomes 
vanishingly small as measured by the TM's mean free path 
A . But, nevertheless, it remains on order of the effective 
"volume of collision" , ~ 7rrQ A . This observation prompts 
that fc-particle correlations are concentrated just at those 
particular subsets of fc-particle phase space which corre- 
spond to collisions. Therefore effects of the correlations 
completely hold under the limit. The same is said by the 
equality §8§ which also holds out. This becomes obvious if 
we take into account that actually important parameter of 



the integrated CF under the Boltzmann-Grad limit must 
be A instead of n and rewrite (JSJ> in the form 

oo k 

W (t, A; A) + ]T ^ W k (t, A; A) = W (t, A; A/(l + «)) , 

k=l 

where W k {t, A; A) = lim n k V k (t, A; n) . Thus, in essence 
nothing changes under the Boltzmann-Grad limit. 

It is necessary to underline that the correlations under 
our attention are qualitatively different from correlations 
which appear in standard approximations of the BBGKY 
hierarchy and connect velocities of particles after colli- 
sion (see e.g. [6]). In our notations, a pair correlation 
of such the kind would look nearly as Vi(t, R, r, P, p) = 
Fi(t, R, r', P', p') - Fi(t, R, r, P, p) , where the primed 
variables describe the pre-collision state. Clearly, because 
of the phase volume conservation during collision, integra- 
tion of this expression over p — r — R and the momenta 
results in zero. This observation shows that correlations 
under our interest do live not so much in momentum space 
as in the configuration space. Their salt is that they con- 
nect coordinates and walks of particles and may coexist 
with statistical independence of particles' velocities. By 
their very nature, they are attributes of spatially inho- 
mogeneous states and evolutions (evolution of Vo(t, A; n) 
gives an example). More profoundly, statistical and phys- 
ical meanings of such correlations were discussed already 
in [7] (so for detail see [7]d). 

In view of the aforesaid0 we have to conclude that 
the Boltzmann-Lorentz equation [2,5] and, moreover, the 
Boltzmann equation in itself do not represent a (low- 
density) limit of the exact statistical mechanical theory. 
The classical kinetics is only its approximate probabilistic 
model (may be good in one respects but caricature in oth- 
ers). Of course, in the exact theory also molecular chaos 
does prevail. But here it is much more rich, even (and 
first of all) in case of dilute gas, and does not keep within 
naive probabilistic models. 
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